UNCLASSIFIED 


r 


Reproduced 
luf.  Ute 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  ■whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


295  787 


62-ijii 


'J.'iJ  iTYPEfiSOii  J.C  AREA  iULE 


l;,  Ladyshenski;/ 


29  573? 


FTP-Tr-62-im/i+2+4 


UNEDITED  ROUGH  DRAFT  TRANSLATION 


THE  HYPERSONIC  AREA  RULE 


By:  M.  D.  Ladyzhenskiy  (Moscow) 
English  Pages:  10 

Source:  Inzhernernyy  Zhurnal,  Vol.  1, 
No.  1,  1961,  pp.  159-163 

[An  abstract  of  this  article 
appears  in  S./124-62-0-2-3/14] 


SC-1585 


THU  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGI- 
NAL  FOREIGN  TEXT  WITHOUT  ANT  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 
ADVOCATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 
AND  DO  HOT  NECESSARILY  REFLECT  THE  POSITION 
OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  DI¬ 
VISION. 


PREPARED  BY« 

TRANSLATION  SERVICES  BRANCH 
FOREIGN  TECHNOLOGY  DIVISION 
WP-AFB,  OHIO. 


FTD-TT 


Pi>t6  3  jan  1^63 


P 


THE  HYPERSONIC  AREA  RULE 
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M.  D.  Ladyzhenskiy  (Moscow) 


The  hypersonic  area  rule  has  been  formulated  in  a 
work  [Ij  and  in  a  hypothesis  [2]  that  all  the  mass  of 
a  gas  is  concentrated  in  an  infinitely  thin  layer 
contiguous  to  the  shockwave.  According  to  that  rule 
when  thin  non-axially-symmetrical  blunt  bodies  with 
the  same  bluntness  drag  values  and  the  same  laws 
of  transverse  area  change  in  the  direction  of  the 
stream,  are  subjected  to  flow,  the  shock-wave  sur- 
Ihces,  the  laws  of  pressure  change,  and  consequently 
also  the  forces  of  resistance  acting  on  the  body 
coincide,  and  in  this  process  the  surfaces  of  the 
shock  waves  have  axial  symmetry. 


i.  Determining  the  results'  bounds  of  applicability  [1].  As 
an  example  of  applying  the  hypersonic  area  rule  [1]  let  us  construct 
a  body  equivalent  to  a  thin  cone  of  revolution,  i.e,  having  the  same 
bluntness  drag  values  and  the  sams  law  of  cross-sectional  area  change 
in  the  stream  direction  as  a  cone.  The  cross-section  of  the  body 

is  assumed  to  be  elliptical,  the  larger  semiaxis  of  which  equals  the 

radius  of  the  shock  wave  and  the  area  of  which  equals  the  area  of  the 

cross-section  of  the  cone  of  revolution  (Pig.  1).  In  other  words  the 

eccentricity  of  the  ellipse  in  each  section  has  the  largest  possible 
value  compatible  with  the  demand  (condition  3  in  work  [1])  that  the 


/ 
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body  must  not  exceed  the  volume  delimited  by  the  shock  wave. 

As  was  noted  In  the  fork  [1 J  the  area  rule  may  be  combined  with 
the  law  of  similarity  In  the  flow  about  thin  bimt  bodies  [2],  as  a 
result  of  which  dimensionless  values  characterizing  the  flow  are  defined 
3-t  a  fixed  value  of  the  index  of  adiabatic  curve  h  by  two  dimensionless 
parameters:  the  known  parameter  of  similarity  of  flow  about  thin 
blunt  bodies  K  =  M  t  and  the  parameter  Kx  =  tt/2c  S)1//2Lt2  character- 
izing  the  effect  of 'bluntness  and  equal  in  order  of  size  to  the 
square  root  of  the  ratio  of  body  drag  to  bluntness  drag.  Here  T~Sly^2/L 
is  a  small  dimensionless  value  characterizing  the  thickness  of  the 
body;  S#  is  a  certain  representative  area  of  the  body's  cross-section; 

L  is  the  length  of  the  body;  cx,  S  are  coefficient  of  bluntness  drag 
and  the  area  of  the  mid-section  of  the  bluntness  [surface  area  per 
unit  depth]  respectively. 

The  shape  of  the  bluntness  is  non-essential  in  the  assumption 
[1]  that  the  effect  of  bluntness  may  be  replaced  by  the  effect  of  a 
blast  at  the  forward  point  of  the  body  with  an  energy  equal  to  blunt¬ 
ness  drag.  In  view  of  this  the  area  of  bluntness  S  is  introduced 
into  the  expression  for  Kx  instead  of  the  diameter  of  the  bluntness. 

Let  us  suppose  that  the  Mach  number  M  of  the  free  stream  equals 
infinity.  Then  when  h  is  fixed  the  dimensionless  variables  will 
depend  upon  the  single  parameter  Ki . 

Figure  2  shows  in  the  form  of  curves  (for  the  case  of  K  =  1.4) 
the  relationships  of  the  values  k  (ratio  of  the  major  to  the  minor 
semiaxis)  and  (X  -  X0)/X0  (ratio  of  body  drag  without  bluntness 
drag  to  bluntness  drag)  in  the  function  Kx  =  (tt/2cxS)^2L  tan2  a  where 
a  is  half  the  apex  angle  of  the  circular  cone.  The  shape  of  the 
shock  wave  was  defined  from  the  solution  of  the  problem  of  flow 
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around  a  thin  blunted  cone  in  accordance  with  Chernyy  [2].  It  is 
sensible  to  use  the  area  rule  when  x/XQ  >  1.1,  which,  as  follows 
from  Pig.  2,  corresponds  to  Ki  ^  0.1.  At  lesser  values  of  Ki  the 
body  drag  is  practically  determined  by  the  value  of  the  bluntness 
drag.  At  large  values  of  Ki  the  area  rule  becomes  invalid  when  k 
approaches  unity  (more  accurately  when  k-1  ~  (h-1)/(h+1)  [1]),  which 
occurs  approximately  at  K*  =  1.2. 

Thus,  the  range  of  application  of  the  area  rule  lies  within  the 
limits  0.1  <  Ki  <  1.2.  Moreover  the  ellipse  in  the  cross-section  of 
the  body  may  have  an  amply  extended  shape  differing  from  a  circle 
(15 k  1.5).  This  result  may  be  of  practical  interest.  But  seeing 
that  the  results  are  obtained  in  rough  assumptions  about  the  concen¬ 
tration  of  all  the  gas's  mass  in  an  infinitely  thin  layer  behind  the 
shock  wave  we  must  define  more  exactly  the  area  rule. 

2.  A  more  exact  definition  of  the  area  rule.  As  before  we 

assume  that  M  >  1,  but  in  contrast  to  another  work  [1],  we  do  not 
00 

impose  the  compulsory  condition  M^T  >1.  We  introduce  a  cylindrical 
system  of  coordinates  xL,  yL,  0  (the  x  axis  passes  through  the  leading 
point  of  the  body  and  is  directed  along  the  stream).  Let  us  designate 
by  uU^,  vU^,  wUm  the  components  of  the  velocity  in  the  axial,  radial, 
and  circular  directions  respectively;  PP^U^  is  the  pressure,  pp^ 
the  density,  U  the  velocity  of  the  free  stream  along  the  x  axis, 
p  its  density.  We  will  characterize  the  bluntness  value  with  the 
dimensionless  bluntness  diameter  dL,  where  d  is  a  small  quantity.  We 
will  write  the  equation  of  the  body  surface  in  the  form  y  =  Tf(x,  0). 

Let  us  separate  the  entropic  layer,  lie.,  the  region  occupied  by 
the  lines  of  flow  which  have  passed  that  area  of  the  surface  of  the 
shock  wave  where  the  angles  of  inclination  formed  by  the  surface  of  the 
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shock  wave  with  the  direction  of  the  free  flow  are  not  small  (Fig.  3). 
Let  the  equation  of  the  provisorily  interpolated  boundary  of  the 
entropic  layer  be  y  =  5cp(x,  0)  where  6  is  a  small  quantity.  Beginning 
at  a  certain  x  =  xQ  ~  d  the  angles  of  inclination  of  the  boundary  of 
the  entropic  layer  will  have  order  6 . 

Below  are  listed  the  estimates  of  the  parameters  of  the  stream  in 
the  entropic  layer  similar  to  those  outlined  in  another  work  J •  Let 
us  also  note  that  the  influence  of  the  entropic  layer  on  the  pressure 
distribution  on  a  thin  blunt  cone  is  also  examined  in  a  work  [4], 

Let  us  assume  that  on  the  surface  of  the  entropic  layer  the 
relationship  p  ~  da  holds  true,  where  a  is  a  positive  number  which 
can  be  determined. 

As  results  from  the  following,  the  order  of  pressure  across  the 

a 

entropic  layer  is  sustained,  therefore  we  may  write  p  ~  dK  for  the 
density,  using  the  conditions  of  adiabaticity .  Let  us  now  write  the 
continuity  equation  for  the  entropic  layer.  Equating  the  flow  rate  in 
the  entropic  layer  with  the  flow  rate  in  the  stream  channel  in  the 
free  stream  across  an  area  equal  to  the  midsection  of  the  bluntness 
we  have 

<!"•  ~  p  u  5,  (  i  ) 

where  a  is  the  area  occupied  by  the  entropic  layer  in  the  cross- 
section  where  x  =  const  (hatched  in  Fig.  J>) .  Insofar  as  this  follows 
from  the  Bernouilli  equation  in  the  entropic  layer  u  ~  1  we  have 
a  =  d2  =  —  .  It  is  obvious  that  for  6  entering  the  boundary  equation 

14/ 

of  the  entropic  layer  we  get  the  evaluation 

6 2 .  .v  o,  (2) 

where  S  is  the  area  of  a  cross-section  of  the  body.  Let  us  require 
that  the  area  of  the  body  in  its  order  of  magnitude  not  exceed  the 
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entropic  layer  S  <  a.  Then,  obviously 


Since  the  usual  estimate  p  ~  6 2  for  the  hypersonic  stream  on  the 
exterior  boundary  of  the  entropic  layer  is  correct  for  pressure  we 
derive  the  equation  for  determining  a: 

<1'  ~  ~  T'  a  ,_2*_ 

x  !-  1 

Finally  for  the  stream  parameters  in  the  entropic  layer  we  have: 

V  ~  rf*7'.  t<  ~  d7*-"1,  6  ~  * *  (  ^  ) 

on  condition  that  t  and  d  are  connected  in  order  of  values  by  the 
relationship  derived  from  the  condition  S  ~  a: 

x 

r  ~  oK^ 

This  relationship  practically  coincides  with  the  condition  t  ~ Vd  > 
expressing  the  fact  that  body  drag  is  comparable  in  order  of  magnitude 
with  bluntness  drag  [2], 

Let  us  estimate  the  pressure  drop  in  the  entropic  layer. 

We  have  from  the  equation  of  motion 

'  f/L  ,  ,  _  0  (U  !“L  +V±+JL  £*), 

0!f  \  tty  y  </0  (*7  ) 

J-  It  ,,  ... (u  !>“:  A.  r  :V.  +  11  '!0L\  . 
y  oi)  \  <u  oy  ;/  01}  J 

Consequently  -  y  ,  since  all  terms  on  the  right  side  of  Eq.  7 
are  of  the  same  order  of  smallness.  Hence,  taking  Eqs.  5  and  6  into 
consideration,  the  following  estimate  is  correct  for  the  drop  in 
pressure  in  the  radial  as  well  as  in  the  circular  direction: 

2(*  :-I> 

Ap  ~~  U"  T  *  . 

Thus  the  pressure  in  the  entropic  layer  may  be  considered  constant 


with  a  relative  error 


-  rf*Tl  _  T  X 


somewhat  greater  than  the  relative  error  in  the  theory  of  slight 


-5- 


disturbances  in  the  hypersonic  stream,  equalling  t2  as  is  well  known. 

Let  us  formulate  the  flow  problem.  When  x  <  xQ  let  there  be  an 

axisymmetrical  nose-piece  of  a  body  with  an  axis  of  symmetry  along 

the  x  axis,  the  flow  about  which  is  completely  calculated.  As  follows 

from  what  has  been  set  forth,  since  the  entropic  layer  cannot  restrain 

the  drop  in  pressure  in  the  circular  direction  the  pressure  in  the 

cross-section  where  x  =  const  in  the  area  where  x  >  x  on  the  exterior 

o 

boundary  of  the  layer  must  be  constant.  To  fulfill  this  condition  it 
is  enough  that  the  surface  delimiting  the  entropic  layer  have  axial 
symmetry.  (its  equation  may  moreover  be  written  y  =  6Y(x).) 

Then  the  flow  outside  of  the  entropic  layer,  axisymmetrical 
according  to  agreement  when  x  <  x  ,  preserves  axial  symmetry  also 
when  x  >  x,  and,  consequently  the  stipulation  that  the  pressure  be 
constant  in  the  circular  direction  on  the  exterior  boundary  of  the 
entropic  layer  will  be  fulfilled. 

Let  us  derive  the  equation  connecting  S  and  o  when  x  >  xQ,  for 
which,  as  in  Sychev  [j5],  we  use  the  continuity  equation.  We  will 
designate  values  in  the  surface  xQ  with  the  subscript  0.  Separating 
the  elementary  stream  channel  in  the  entropy  layer  we  will  write  this 
flow-rate  equation  for  it: 


Pii'myurlQod!/',  „  puydQdy. 


(9) 


From  adiabaticity  equations  of  Bernouilli,  too,  who  rejects  terms 
of  order  t2  we  have  for  p  and  u: 


f>  =  f'O 


— >  *  a*. 

-  x  —  i  p0  \  p„  I 


1 

(x  -  1)  At- 


(10) 


Let  us  analyze  Eq.  9  for  pu  and  integrate  over  the  whole  area  of  the 
entropy  layer  in  the  cross-section  x  =  xo  (we  utilize  the  circumstance 
that  the  coordinates  of  the  chosen  current  lines  y  and  0  satisfy  the 
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relationships  y  =  y(y  ,  ®0)>  0  -  0(yo>  0O))*  0n  the  right  side  of  the 
equality  we  obviously  will  obtain  area  a  occupied  by  the  entropy  layer 
in  cross-section  x*  Taking  into  consideration  the  axisymmetry  of  the 
boundary  of  the  entropy  layer  we  have  a  =  7r02Y2(x)  -  S.  Finally,  the 
desired  relationship  is  written 


=  (r)  -  F  {;.)  ..  S  W. 


F(P) 


-  ss  m 


1  -J- _ ” _ L_  _  Pt_ 

*  —  1  .1/^  X  -  1  Po 

1  +  _JL _ 1 _ Pt,  /£\~ 

*-l  ,1/*  x  — 1  Po  \Pa  J 


.  yjibtplyo- 
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Axisyxnmetrical  flow  outside  the  mtropy  layer  when  x  >  xQ  may  be 
computed  by  one  of  the  known  exact  methods,  e.g.,  by  the  method  of 
characteristics.  The  boundary  of  the  entropy  layer  will  be  defined 
moreover  in  the  process  of  solving  from  Eq.  11  which  plays  the  role 
of  an  edge  condition  replacing  the  condition  of  non-flow.  Hence  it 
follows  that  flow  when  x  >  xQ  is  completely  determined  by  the  assign¬ 
ment  of  a  law  of  change  in  the  area  of  the  cross-section  of  the  body 
S(x).  In  this,  one  must  impose  one  more  obvious  limitation  on  the 
shape  of  the  body:  it  cannot  exceed  the  limits  of  the  "entropy  circle" 
(see  Fig.  3),  which  may  be  symbolically  written 

p'  c  (x).  (12) 

We  may  now  formulate  the  more  precise  hypersonic  area  rule  'in  the 
following  way.  When  thin  blunt  bodies  having  axisymmetrical  nose- 
pieces  coinciding  for  some  distance  from  the  leading  point  of  the 
body  and  having  also  the  same  laws  of  change  in  the  areas  of  the 
cross-section  of  the  remaining  parts: 

a)  flows  outside  of  the  entropy  layers  are  axisymmetrical;  the 
parameters  of  the  flows  at  corresponding  points,  the  surfaces  of  the 
shock  waves,  and  the  provisorily  introduced  boundaries  of  the  entropy 
layers  coincide; 
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b)  pressure  in  the  entropy  layers  depends  only  on  x  and  the  law 
of  pressure  change  is  the  same  for  the  bodies  under  investigation; 
as  a  result  of  this  the  resistance  forces  acting  on  the  bodies  are 
equal  as  the  resistance  X  is  expressed  in  the  form 

-V  =  -V.  +  $  V  U)  ,,  <*)  ,U.  (  15  ) 

where  XQ  is  the  resistance  of  the  nosepiece  of  the  body.  In  addition 
it  is  assumed  that  conditions  6  and  12  are  fulfilled. 

The  results  obtained  are  easily  generalized  to  the  case  of  flows 
with  dissociation.  A  consideration  of  these  phenomena  leads  only  to 
a  change  in  the  form  of  function  F(p)  in  Eq.  11. 

5.  Comparison  of  the  results.  Let  us  compare  the  result  we 
have  obtained  with  the  area  rule  proven  in  another  work  [lj.  The 
requirement  for  the  coincidence  of  the  laws  of  change  in  area  of  a 
cross-section  in  the  direction  of  the  x  axis  and  also  the  stipulation 
that  body  drag  must  not  exceed  bluntness  drag  in  order  of  magnitude  are 
common  to  both  theorems.  The  difference  in  the  theorem  proven  in  section 
2  in  the  formulation  consists  in  the  following: 

a)  the  obligatory  condition  Mw  t  >  1  is  not  imposed  as  is  done 
in  the  other  work  [1]; 

b)  Instead  of  the  requirement  for  equality  in  the  magnitudes  of 
bluntness  resistance  [1]  a  more  severe  limitation  is  imposed:  the 
nosepieces  of  the  equivalent  bodies  must  coincide  in  form,  being 
axisymmetrical; 

c)  instead  of  condition  5  in  the  work  [1], according  to  which  the 
body  does  not  exceed  the  volume  delimited  by  the  surface  of  the  shock 
wave, the  severer  limitation  12  is  imposed  that  the  body  must  not  exceed 
the  limits  of  the  inner  boundary  of  the  entropy  layer. 
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As  a  result  of  this  we  may  expect  that  the  values  of  k  found  In 
section  1  characterizing  the  difference  of  cross-section  of  the  body 
from  the  cross-section  of  the  equivalent  body  of  revolution  are  some¬ 
what  too  high. 


Pig.  1. 

1) Shock  wavej  2) 
cross  section  of 
the  circular  body; 

5)  cross  section  of 
the  equivalent  body. 


Fig.  5. 
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